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We have attempted to find the global minima of clusters containing between 20 and 80 atoms 
bound by the Morse potential as a function of the range of the interatomic force. The effect of de- 
creasing the range is to destabilize strained structures, and hence the global minimum changes from 
icosahedral to decahedral to face-centred-cubic as the range is decreased. For A*' > 45 the global 
minima associated with a long-ranged potential have polytetrahedral structures involving defects 
called disclination lines. For the larger clusters the network of disclination lines is disordered and 
the global minimum has an amorphous structure resembling a liquid. The size evolution of polyte- 
trahedral packings enables us to study the development of bulk liquid structure in finite systems. 
As many experiments on the structure of clusters only provide indirect structural information, these 
results will be very useful in aiding the interpretation of experiment. They also provide candidate 
structures for theoretical studies using more specific and computationally expensive descriptions 
of the interatomic interactions. Furthermore, Morse clusters provide a rigorous testing ground for 
global optimization methods. 
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I. INTRODUCTION 



Structural information is of fundamental importance in 
addressing the chemical and physical properties of any 
system. Unfortunately, there is no direct experimen- 
tal method for determining the structure of free clus- 
ters in molecular beams. Instead, one measures prop- 
erties which depend on structure and employs mod- 
els of the predicted favoured geometries. This ap- 
proach has been|-combined with techniques such as elec- 
tron diffraction,EJ mass, spectral abundances^ ^chemical 
reactivityp magnetismcl and x-ray spectroscopyJj The in- 
version of the experimental data to obtain structural in- 
formation, though, can be problematic, and always relies 
on comparisons with the predictions of structural models. 

Sizes exhibiting special stability are known for^certain 
morphologies and classes of interatomic potential.H'Q The 
size-dependence of properties sometimes reveals these 
magic numbers, and thus enables a confident structural 
assignment to be made. However, often rather little is 
known about the structure between these magic num- 
bers. One of the most powerful experimental techniques 
that addresses this deficiency is the flow-reactor approach 
which probes the chemical reactivity of size-selected clus- 
ters. For example, this method has been applied to 
nickel clusters, to. give detailed information for all sizes 
up to Af^Tl.rB The results show that around A^=13 
and A^=55 the clusters are icosahedral, in agreemen|t.jwjtb, 
the observed magic numbers in other experiments.aii3't3 
However, in the size range 29 < TV < 48 only one struc- 



tural assignment has so far been made because of the 
large number of possible geometries to be considered and 
the presence of multiple isomers. til 

The theoretician can aid in the task of structural as- 
signment by providing realistic candidate structures. In- 
deed, many studies have modelling specific clusters, but 
ah initio calculations are only feasible for small sizes, es- 
pecially for transition metals, and so empirical potentials 
are often used. However, as is clear from the dissepity of 
theoretical results obtained for nickel clusters JlsEB con- 
sensus between methods is lacking and it is hard to know 
which (if any) of the results should be believed. Even 
with the simplified description of the interatomic inter- 
actions provided by an empirical potential, it can be an 
extremely difficult task to search the potential energy 
surface (PES) extensively enough to be confident that 
the global minimum has been found. Also, many empir- 
ical potentials are too complicated to provide an under- 
standing of the relationship between the potential and 
the observed structure and so little physical insight is 
gained. 

Therefore, to understand cluster structure there is a 
need for a hierarchy of theoretical models from the gen- 
eral to the specific. In the present study we use a simple 
model to understand the structural effects of the range of 
the potential, and so provide one part of the framework 
for understanding the physical basis of cluster structure. 
We are confident that we have found most of the global 
minima giving the most comprehensive model of cluster 
structure in the small size regime. 
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FIG. 1. (a) 38-atom truncated octahedron, (b) 55-atom 
Mackay icosahedron, and (c) 75-atom Marks decahedron. 
These clusters have optimal shapes for the three main types 
of ordered packing seen in clusters: face-centred cubic (fee), 
icosahedral and decahedral, respectively. The latter two 
morphologies cannot be extended to the bulk because of the 
hve-fold axes of symmetry. 



One of the most interesting aspects of cluster structure 
is the manifestation of non-crystallographic symmetries 
which arise from the absence of translational periodicity. 
In particular, many clusters are found to have fivefold 
axes of symmetry, including two of the three main types 
of ordered structure adopted by simple atomic clusters. 
Decahedra have a single fivefold axis of symmetry and 
are based on pentagonal bipyramids, while icosahedra 
have six fivefold axes of symmetry. The third morphol- 
ogy consists of close-packed clusters. Particularly stable 
examples of each type are illustrated in Fig. ^. 

All the above morpholgies have been observed exper- 
imentally. Many gas phase clusters have been shown to 
be icosahedral through the presence of the noagic num- 
bers associated witiL-th^ Mackajj— imsahedraO in mass 
spectra: j-caxe gasescHtJ, metals, oEJlZI and molecular 
clusters. ESc^ Icosahedral and decahedral structures are 
also conuaonly reported for metal clusters supported on 
surfaces,E3 and more recently fee and decahedral clus- 
ters have been_pbserved for gold clusters passivated by 
alkylthiolates.El"El 

All three structural types are also exhibited by 
Lennard-Jones (LJ) clusters. For N < 1600 icosahedra 
are most stable; from this size up to N Ki 10L.decahedra 
are most stable and above this fee clusters.E3 However, 
these changes do not occur abruptly. The glpbaL mini- 
mum of LJ38 is the fee truncated octahedroi£3S (Fig. 
||(a)) and for at least six sizes with < 110 tijc_£lobal 
minimum is based upon a Marks deeahedroioLj (Fig. 
|(c)). In this paper we consider a potential with variable 
range to provide a model system which exhibits a much 
greater diversity of structural behaviour than LJ clusters 
in the small size regime. Consequently, the results are 
relevant to a much wider range of systems. 

There have been a number of previous studies on the 
effect of the range of the potential on. the structure and 
phase behaviour of small elusters.EaeJ These have shown 
that the number of minima and saddle points on the PES 
increases ja^i-tlie, range decreases — the PES becomes more 
ruggedEj'OcElo — and strained structures are destabi- 
lized. The latter effect results in range-induced transi- 
tions between the ordered morphologies ,E3 and the desta- 



lance of the liquid phase as the 
3 In a previous paper we exam- 



bilization and disaxi 
range is decreased. E 
ined Morse clusters caataining up to 25 atoms and a se- 
lection of larger sizes.E3 Here we consider Morse clusters 
in the size range 20 < iV < 80 atoms. Some of the lowest 
energy structures given here supersede the results of the 
previous paper. 



II. METHODS 



A. The potential 



The Morse potential!^ may be written as 

Vm = e^e''«'i-''"^/''°)(e''°(i-'''^/''°) - 2), (1) 



where e is the pair well depth and tq the equilibrium 
pair separation. We denote an TV-atom cluster bound 
by the Morse potential as M^r. In reduced units (e = 1 
and rp = 1) there is a single adjustable parameter, poj 
which determines the range of the interparticle forces. 
Fig. H shows that decreasing po increases the range of 
the attractive part of the potential and softens the re- 
pulsive wall, thus widening the potential well. Values of 
Po appropriate to-a, range of materials have been cata- 
logued elsewhercElJ The LJ potential has the same cur- 
vature at the bottom of the well as the Morse potential 
when Po = 6. Girifalco has obtained an intermolecu- 
lar potential for Cgo moleculesEll which is isotropic and 
short-ranged relative to the equilibrmrn pair separation, 
with an effective value of po = 13.62E3 The alkali metals 



have longer-ranged interactions, for example po=3.15 has 
been suggested for sodium.lE-Fi 
a value of po=3.96 for nickel.Eil 



'"itting to bulk data gives 
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FIG. 2. The Morse potential for different values of the 
range parameter po as indicated. 
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At absolute zero the structure with the lowest free en- 
ergy is simply the global potential energy minimum of the 
PES. At higher temperatures, entropic factors must also 
be considered. Although we only perform a comprehen- 
sive survey of the zero Kelvin geometries, the structural 
effects of temperature are considered in the subsequent 
discussion. 

To understand the structural effects of the range pa- 
rameter, it is instructive to look more closely at the 
form of the potential. The energy can be partitioned into 
three contributions: 

Vm — —nnn£ + -E'strain + Ennn- (2) 

The number of nearest-neighbour contacts, n„„, the 
strain energy, £^strain, and the contribution to the energy 
from non-nearest neighbours, are given by 

i<j,Xij<Xo 

i^strain = £ (e""""" " 1)"' 

i<j;Xij<xa 

Ennn = e ^ e""""- (e-"""- - 2) , 

i<j,Xij>XQ 



and are approximately spherical, the icosahedra have 
the largest n„„. Complete Mackay icosahedra occur at 
N — 13, 55, 147, ... A pentagonal bipyramid has only 
{111} faces, but because it is not very spherical more 
stable decahedral forms are obtained by truncating the 
structure parallel to the five- fold axis to reveal five {100} 
faces and then introducing re-entrant {111} faces be- 
tween adjacent {100} faces. The resulting structure is 
called a Marks decahedron (Fig. ||(c)) and was fiEsi. pre- 
dicted by the use of a modified Wulff construction. l£| Dec- 
ahedra generally have lower values of n„„ than icosahe- 
dra because of the {100} faces. The tetrahedron and 
octahedron are fee structures that have only {111} faces, 
but they are not very spherical. The optimal fee struc- 
ture is the truncated octahedron with regular hexago- 
nal faces (Fig. ||(a)). At larger sizes than those consid- 
ered in this study, the optimal structure involves further 
facetting, SQ_Lhat it more closely approximates the Wulff 
polyhedron.Ell Of the three morphologies fee structures 
have the smallest values of n„„. 



(3) 




where 



i/fQ ~ 1, and xq is a nearest- neighbour 



criterion. Xij is the strain in the contact between atoms 
i and j. 

The dominant term in the energy comes from n„„. 
Ennn is a Smaller term and its value varies in a sim- 
ilar manner to n„„. It is only likely to be important 
in determining the lowest energy structures when other 
factors are equal. For example, bulk fee and hexago- 
nal close-packed (hep) lattices both have twelve nearest 
neighbours per atom. Next-nearest neighbour interac- 
tions are the cause of the lower energy of the hep ccyglal 
when a pair potential such as the LJ form is usedO'Lj 

^-strain, which measurcs the energetic penalty for the 
deviation of a nearest-neighbour distance from the equi- 
librium pair distance, is a key quantity in our analysis. 
It must not be confused with strain due to an applied ex- 
ternal force. For a given geometry, i^strain grows rapidly 
with increasing po because the potential well narrows. To 
a first approximation the strain energy grows quadrati- 
cally with po-t3 Hence, decreasing the range destabilizes 
strained structures. 

From the above analysis we can see that minimiza- 
tion of the potential energy involves a balance between 
maximizing 7i„„ and minimizing i?strain- The interior 
atoms of the three morphologies (Fig. |) all have twelve 
nearest neighbours, and so differences in n„„ arise from 
surface effects. The optimal shape for each morphology 
results from the balance between maximizing the pro- 
portion of {111} faces (an atom in a {111} face is 9- 
coordinate, but in a {100} face only 8-eoordinate) and 
minimizing the fraction of atoms in the surface layer. 
As Mackay icosahedra (Fig. |(b)) have only {111} faces 





FIG. 3. Examples of the strain involved in packing tetra- 
hedra. (a) Five regular tetrahedra around a common edge 
produce a gap of 7.36°. (b) Twenty regular tetrahedra about 
a common vertex produce gaps equivalent to 1.54 steradians. 



A Mackay icosahedron can be decompos ed in to twenty 
fee tetrahedra, 



but as we see from Fig. II A (b), when 
twenty regular tetrahedra are packed around a common 
vertex large gaps remain. To bridge these gaps the tetra- 
hedra have to be distorted, thus introducing strain. The 
distance between adjacent vertices of the icosahedron 
is 5% larger than the distance between a vertex and 
the centre. Similarly, a pentagonal bipyramid can be 
decomposed into five fee tetrahedra sharing a common 
edge. Again, a gap remains if regular tetrahedra are 



used (Fig. II A (a)) and consequently decahedral struc- 



tures are strained, although not as much as icosahedra. 
In contrast, close-packed structures can be unstrained. 

Having deduced the relative values of n„„ and i?strain 
for icosahedral, decahedral and fee structures, we can 
predict the effect of the range on the competition between 
them. For a moderately long-ranged potential, the strain 
associated with the icosahedron can be accommodated 
without too large an energetic penalty and so such struc- 
tures are the most stable. As the range of the potential is 
decreased, the strain energy associated with icosahedra 
increases rapidly, and there comes a point where deca- 
hedra become more stable. Similarly, for a sufficiently 
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short-ranged potential fee struetures beeome more stable 
than deeahedral struetures. 

The above decomposition of the potential energy also 
helps us to understand the effect of size on the energetic 
competition between the three morphologies. The dif- 
ferences in n„„, which arise from the different surface 
structures, are approximately proportional to the sur- 
face area (oc N^^^). The strain energies, however, are 
proportional to the volume (oc N). Therefore, the dif- 
ferences in ii^strain increase more rapidly with size than 
the differences in n„„, thus explaining the change in the 
most stable morphology from icosahedral to decahedral 
to fee with increasing size. The effect of increasing the 
size is similar to the effect of decreasing the range of the 
potential: both destabilize strained structures. 



B. Searching the potential energy surface 

The principal method that we have used to generate 
candidate structures for the global minima makes use 
of the physical insight gained from the last section. We 
have simply attempted to construct geometries that max- 
imize n„„ for the three ordered morphologies .0 The re- 
sulting structures were then optimized hjj either conju- 
gate gradienttfl or eigenvector-foUowingEj methods. A 
similar approach was successfully used by Ncwthby to 
generate low energy icosahedra for LJ clusters .eI The ef- 
fectiveness of this method is demonstrated by how few of 
Northby's lowest energy structures have been superseded 
and by theJpHffitkwpf-tjme that it has taken to find these 
exceptions .QO^E3'E2rE3 Our approach, however, depends 
on the imagination of the practitioner to conceive of all 
the possible ways that a structure with a large value of 
n„„ could be obtained. Furthermore, this method will 
always fail to find the global minimum if the latter is not 
based on one of the ordered structures, as is the case for 
the larger clusters we have considered at low values of po- 

To complement the above approach two global opti- 
mization techniques were used to try to find structures 
that might have been missed. Firstly, we used a method 
in which eigenvector-following is emn lovR d to take steps 
directly between minima on the PES.E3H If low temper- 
ature Metropolis Monte Carlo is used in this space of 
minima, the system will walk down to the bottom of a 
basin containing many minima. This technique avoids 
the difficulties associated with trapping in local minima 
that can occur for methods which take steps directly in 
configuration space. Secondly, we used a 'basin hop- 
ping' approach, which has proved to be effective for LJ 
clusters;E3E3 it is the only unbiased global optimization 
method to have found the global minima that are Marks 
decahedra. For Morse clusters it was able to find all the 
lowest energy minima at /Oo=3, 6, 10 and all but twelve 
at po=14; this included some structures that were lower 
than any of those we had constructed. The fact that we 
found most of the minima both by unbiased global op- 



timization and by construction makes us confident that 
our lowest energy structures are truly global minima. 

The results for the basin-hopping algorithm are im- 
pressive because the global optimization task for Morse 
clusters is a difficult one. The size of the configuration 
space for the larger clusters is compounded for short- 
ranged potentials by the nature of the PES. Firstly, the 
PES becomes more rugged as the range is decreased.CJ 
The physical reason for the larger number of minima at 
short range is the loss of accessible configuration space 
as the potential wells become narrower, producing bar- 
riers whecfi there are none at long range. Stillinger and 
Stillingeil£3 and Bytheway and Kepertc2l both found that 
minimizations performed from random starting configu- 
rations are much less likely to find the global minimum 
for a short-ranged potential. 

Similarly, barrier heights are also likely to become 
higher and rearrangements more localized for a shorter- 
ranged potential. These trends have been observed in 
comparisonajjf the rearrangements of 55-particle Cgo and 
LJ clusters.E3 Hence the range of the potential is likely 
to have a significant impact on the dynamics, making es- 
cape from local minima much more difficult. This effect 
has been observed by Rose and Berry who have shown 
that the rate at which the ground state structure of a 
potassium chloride cluster is found upon cooling can be 
significantly decreased by using a shielded Coulomb po- 
tential to reduce the range of the interactions .EJ 

Furthermore, as a result of the competition between 
the decahedral and the various types of close-packed 
structures for short-ranged potentials, it is more likely 
that there are a number of low energy minima which are 
very close in energy but are structurally dissimilar. Each 
of these minima lie at the bottom of their own funnel 
on the PES. This multiple funnel topography can lead to 
cases where optimization is extremely difficult because 
the free energy barriers for transitions between the fun- 
nels can be large, thus leading to trapping. The worst 
cases are when relaxation down the PES preferentially 
takes the system intp_a funnel which does not end at 
the global minimumjij and when the global minimum 
only becomes thCj-state with the lowest free energy at 
low temperatures.E3 

Finally, fee and decahedral minima are more struc- 
turally dissimilar from minima typical of the liquid-like 
state than the icosahedral structures. Therefore, the 
paths from the liquid to these structures are likely to 
be fewer and loiiger than those leading to the icosahe- 
dral structures,E3 and so relaxation down the PES from 
the liquid-like state to fee and decahedral global minima 
is harder. 

For a number of reasons Morse clusters are an ideal 
system with which to test global optimization methods 
designed for configurational problems. Firstly, in this pa- 
per we provide very good estimates for the energies of the 
global minima. Secondly, this system represents a much 
more general — the global minima have a variety of struc- 
tural types — and tougher examination than is provided 
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by LJ clusters, a muciuused test system for global op- basin-hopping algorithm provide a benchmark that any 
timization algorithms.l£3 And finally, the results for our would-be global optimization method should try to beat. 
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FIG. 4. Zero temperature 'phase diagram' showing the variation of the lowest energy structure with A'^ and po- The data 
points are the values of po at which the global minimum changes. The lines joining the data points divide the phase diagram 
into regions where the global minima have similar structures. The solid lines denote the boundaries between the four main 
structural types — icosahedral, decahedral, close-packed and those associated with low po (L) — and the dashed lines are internal 
boundaries within a structural type, e.g. between icosahedra with Mackay and anti-Mackay overlayers, or between decahedra 
with different length decahedral axes. 



III. RESULTS 

All the global minima that we have found are cata- 
logued in Table I along with their energies, point groups, 
number of nearest neighbours, strain energies and the val- 
ues of Po for which they are probably the lowest energy 
minimum. 

The results are summarized in Fig. ^ which provides 
a zero temperature 'phase diagram', showing how the 
global minimum depends upon TV and po- The struc- 
tural behaviour of Morse clusters with fewer than eight 
atoms is rather uninteresting because the global mini- 
mum is independent of po. For all N > 8, however, the 
global minimum changes at least once as a function of 



Pq. For > 13, icosahedral, decahedral and fee struc- 
tures all exist, and the form of the phase diagram is 
in good agreement with the predictions we made earlier 
based upon the decomposition of the potential energy. 
For most sizes the structure changes from icosahedral to 
decahedral to close-packed as the range of the potential 
is decreased. For N < 23, however, the transition from 
decahedral to close-packed occurs at larger values of po 
than we consider in this study. There are also a number 
of sizes (A=38-40, 52, 53, 59 and 61) for which there is a 
transition directly from an icosahedral to a close-packed 
structure; this occurs when ?i„„ for the lowest energy 
close-packed structure is greater than or equal to that 
for the lowest energy decahedron. 
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FIG. 5. Plots of A2-E; as a function of N for (a) po^3, (b) 
po=6, (c) po=10 and (d) po=14. 



The boundaries between the different morphologies are 
sensitive functions of N. Such size dependence is ob- 
served for many properties of clusters, and gradually 
lessens as N increases (because the addition of a single 
atom becomes a smaller perturbation) until the bulk limit 
is reached. The decahedral to close-packed boundary is 
particularly sensitive, because the range of po for which 
the decahedron is most stable changes dramatically even 
when the difference in n„„ between the decahedral and 
close-packed structures changes by one. 

Sizes for which a morphology is the lowest in en- 
ergy for a particularly large range of po indicate that 
the structure is especially stable. The optimal geome- 
tries shown in Fig. || are good examples. Another in- 
dicator of special stability is provided by A2E{N) — 
E{N + + E{N - 1) - 2E{N). Peaks in Aa^; have been 
found to correlate well with the magic numbers (sizes 
at which clusters are particularly abundant) observed in 



Iff for a 



mass spectra.O Plots of /S.2E are shown in Fig. 
number of values of po- Unsurprisingly, the plot for po—Q 
is very similar to that for LJ clusters with peaks due to 
especially stable icosahedral clusters. At higher values 



of po peaks corresponding to close-packed and decahe- 
dral clusters begin to occur. The plot at pn—U is very 
similar to that recently obtained for Ccq clusters using 
the Girifalco intermolecular potentialO If the energy is 
'normalized' by subtracting a suitable function of N, par- 
ticularly stable sizes again stand out (Fig. III). 
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FIG. 6. Comparison of the energies of icosahe- 
dral (solid line with diamonds), decahedral (dashed 
line with crosses) and close-packed (dotted line with 
squares) Mn clusters at po=6. The energy zero is 
Em I, the interpolated energy of Mackay icosahedra. 
Emi = -3.0354 + 0.2624iV^''^ + 8.8400iV^''^ - 6.8381iV and 
was obtained by fitting to the first four Mackay icosahedra 
(iV=13, 55, 147 and 309). 



In the following subsections we will look at the growth 
sequences for each morphology in more detail. We also 
examine the unusual structures that occur for the larger 
clusters at low po, which, as we will see, involve a mixture 
of order and disorder. 



A. Icosahedral clusters 

Many small clusters are polytetrahedral in the sense 
that the whole of the cluster can be divided into tetra- 
hedra. This category includes the 13-atom icosahedron, 
which can be decomposed into twenty tetrahedra shar- 
ing a common vertex. Addition of atoms to the icosa- 
hedron can occur in two ways and the t wo ty pes of 



overlayer that result are illustrated in Fig. Ill A, One 



growth mode (fcc-like) continues the fee packing of the 
twenty strained tetrahedra making up the icosahedron, 
and leads to the 55-atom Mackay icosahedron (Fig. |(b)). 
This scheme introduces octahedral interstices, and so 
the resulting structures are no longer polytetrahedral. 
The other 'anti-Mackay' (hep-like) growth mode involves 
sites which are hep with respect to the tetrahedra. For 
growth on the 13-atom icosahedron, this overlayer pre- 
serves polytetrahedral character. Each of the vertices 
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of the original icosahedron becomes icosahedrally coordi- 
nated, and the structure that results from the completion 
of this overlayer, 45A, is a rhombic tricontahedron; it is 
an icosahedron of interpenetrating icosahedra. Interest- 
ingly, the rhombic tricontahedron is the face-dual of the 
truncated icosahedron made famous by Cgo; indeed, it is 
even a particularly stable shell for the coverage of a Ceo 
molecule by alkaline earth metal atomsEj In previous 
studies, the anti-Maf;kay overlayer has been referred to as 
the polyicosahedralEl or the face-capping overlayer.! Such 
names are reasonable for growth on the 13-atom icosa- 
hedron, but are confusing for growth on larger Mackay 
icosahedra. 



Anti-Mackay 



Mackay 




Second Shell 





Third Shell 




FIG. 7. Atomic positions for the two possible overlayers 
of the icosahedron, anti-Mackay (left) and Mackay (right). 
These are shown for a single face of the icosahedron. 

The icosahedral structu res wi th an anti-Mackay over- 
layer are illustrated in Fig. Ill A and those with a Mackay 
overlayer in Fig. ^ Growth from the 13-atom icosahe- 
dron begins in the anti-Mackay positions, because these 
do not include any low-coordinate edge sites, thus giving 
a larger n„„. However, at some critical size the Mackay 
overlayer becomes lower in energy because of the larger 
strain energies associated with the anti-Mackay overlayer. 
Further growth then leads to the next Mackay icosahe- 
dron. The size at which this change occurs depends on 
the range of the potential; it is given by the dashed line 
of negative slope that divides the icosahedral region of 
the phase diagram (Fig. ^). The crossover size increases 
with the range of the potential. 

At po = 6 the Mackay overlayer is the lowest in energy 
for > 32; the corresponding result for LJ clusters is 
A^ > 31. In a6 initio molecular dynamics calculations for 
lithium clustersjZJ the anti-Mackay overlayer is lowest in 
energy up to A^ = 45. Similarly, polytetrahedral struc- 
tures were observed in a tight-binding study of sodium 



clusters up to the largest size considered A'^=34 (for which 
34 A was the most stable structure) JI3 These effects are 
nicely explained by the long range of alkali metal poten- 
tials. 



21A 21B 



22A 23A 




FIG. 8. Icosahedral global minima formed by growth of 
an anti-Mackay overlayer on the 13-atom icosahedron, ex- 
cept structure 23B which is composed of two face-sharing 
icosahedra. 



Especially stable structures with an anti-Mackay over- 
layer occur when the icosahedral coordination of a vertex 
atom is complete. These structures give rise to the peaks 
in A2-B at A^=19, 23 26, 29, 34 and 45 for po=3; only 
the first four of these peaks are seen at po=6, and only 
the first at po=10 (Fig. |l|). Some of these magic num- 
bers rhf-yEi been observed in—tiic mass spectra of noble 
gase£3~t3 and even bariumJij The centres of the icosa- 
hedra in these structures form a dimer for 19A, an equi- 
lateral triangle for 23A, a tetrahedron for 26A, a trigonal 
bipyramid for 29A, a pentagonal bipyramid for 34A, and 
an icosahedron for 45A. 

At a number of sizes there is more than one global min- 
imum with an anti-Mackay overlayer for different values 
of jOq. The transitions between these structures are re- 
lated to small differences in the value of -Bstrain- Also 



illustrated in Fig. Ill A is structure 23B, which is made 



up of two face-sharing icosahedra. It can be formed from 
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structure 17B by the addition of six atoms to part of the 
over layer. 

The first global minimum with a Mackay overlayer oc- 
curs at N—27. Especially stable structures occur at those 
sizes for which complete faces of the 55-atom Mackay 
icosahedra are missing. These structures give rise to the 
peaks at A^=39, 46 and 49 for po—6 and correspond to 5, 
2 and 1 missing faces, respectively. Again, these magic 
numbeiB have been seen in the mass spectra of noble 



gases.! 




peaks at iV=58, 61, and 64, which correspond to com- 
plete coverage of one, two and three faces, respectively. 

Structure 69C has a vertex atom missing from the un- 
derlying Mackay icosahedron like 38A. There are also 
structures (62C, 65D, 72C and 75B) where an atom is 
added to the surface of the overlayer rather than to the 
Mackay icosahedron. Again there are transitions between 
different anti-Mackay structures resulting from small dif- 
ferences in -^strain- 

The first structure in this size range with a Mackay 
overlayer is 78D. We expect that the crossover from an 
anti-Mackay to a Mackay overlayer will again shift to 
larger size as the range increases, but we have not inves- 
tigated this prediction in the present work. There are 
also three icosahedral global minima, 69B, 70C and 71B, 
which do not fit neatly into either the Mackay or anti- 
Mackay category. Their surface layers have a Mackay-like 
character, but are not in correct alignment with the un- 
derlying icosahedron. The overlayer has been given a 
twist about one of the fivefold axis in order to collapse 
some of the half octahedra at the edge of the overlayer 
into trigonal bipyramidsniji a multiple diamond-square- 
diamond (DSD) process.EZl In fact for M71 at po=5, the 
lowest energy Mackay structure is a C^y transition state 
corresponding to a multiple DSD rearrangement between 
two permutational isomers of 71B. 

Fig. Ill illustrates the variation of the icosahedral en- 
ergies with size. The complete Mackay icosahedra ap- 
pear as narrow minima separated by broad maxima cor- 
responding to structures with incomplete outer shells. At 
Pq—6 it is only near the top of these maxima that fee and 
decahedral structures begin to have similar energies, for 
example at iV=38 and 75. 



FIG. 9. Icosahedral global minima formed by growth of a 
Mackay overlayer on the 13-atom icosahedron. 

Structure 38E has an atom missing from a vertex of 
the original 13-atom icosahedron to allow the overlayer 
to complete a particularly stable form. This structure 
was not found in Northby's stttdy, but it is the lowest 
energy LJ icosahedral cluster .e3 

The icosahedral global minima with more than 55 
atoms are shown in Fig. As for the 13-atom icosahe- 
dron, growth initially occurs at the anti-Mackay sites, be- 
cause this results in structures with larger n„„. Comple- 
tion of this overlayer occurs for a cluster with 127 atoms. 
The vertices of the 55-atom Mackay icosahedron become 
icosahedrally coordinated, but the edge atoms have a 
decahedral coordination shell (this is clearly visible for 
59D) leading to the half octahedra that are visible in 
the surface layers of the anti-Mackay clusters. The most 
prominent peak in A2E due to an anti-Mackay structure 
occurs at 7V=71 for po=6. This corresponds to an over- 
layer which completely covers the five faces surrounding a 
vertex of the underlying icosahedron. There are smaller 



B. Decahedral clusters 
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The decahedral global minima are shown in Figs 
The structures have been grouped according to the num- 
ber of atoms along the fivefold axis of the pentagonal 
bipyramid upon which they are based, and the decahe- 
dral region of the structural phase diagram (Fig. I) has 
also been subdivided on this basis. 

Decahedral clusters grow by capping exposed {100} 
faces and filling in the grooves produced by the re- 
entrant {111} faces. As this process progresses the struc- 
ture changes from prolate to approximately spherical to 
oblate. This cycle begins again when a prolate cluster 
with a longer decahedral axis becomes lower in energy 
than the oblate cluster (e.g. at TV ~ 30 and 54). For 
the clusters based upon a pentagonal bipyramid with 5 
atoms along the fivefold axis {N > 54), the growth pro- 
ceeds asymmetrically — the decahedral axis does not al- 
ways pass through the center of the cluster. For example, 
for 54C the surface structure of the 75-atom Marks dec- 
ahedron is completed on one side of the cluster before 
atoms are added to the other. 
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FIG. 10. Icosahedral global minima formed by growth from the 55-atom Mackay icosahedron. 




FIG. 11. Global minima based upon a decahedron with 
three atoms along the fivefold axis. 



Deviations from this basic growth scheme occur for 
7V=21-30 (Fig. ) and 7V=48-52, 58, 60 and 62 (Fig. [H). 
These structures are formed by addition of atoms to the 
{111} faces surrounding the fivefold axis in sites which 
are hep with respect to the five fee tetrahedra that make 
up the decahedra. These structures are more favourable 
even though they are more strained than the usual dec- 
ahedra because they have a larger n„„. For A^=21-30 
these structures are actually fragments of the 55-atom 
Mackay icosahedron. 

The complete Marks decahedron, 75C, is particularly 
stable. The value of po at which it becomes the global 
minimum (5.81) is the lowest of any of the decahedra. As 
this value-of pmin suggests, it is also the global minimum 
for LJ75.E3 This stability is also indicated by the large 
peak in A2E for /Oo=10 and 14. Other particularly sta- 
ble structures occur at A^=64 and 71; these are fragments 
of 75C with 3 and 4 {100} faces of the Marks decahedron 
complete. 

Decahedral structures have been regularly seen in sup- 
ported metal clusters. Lj However, it is only recently that 
further experimental evidence for the existence of Marks 
decahedra has been founcLiiji [Studies of gold clusters pas- 
sivated by alkylthiolates.Ej"E3 Whetten and coworkers 
were able to isolate fractions with narrow size distribu- 
tions which corresponded to the 75-, 101- and 146-atom 
Marks decahedra. It is significant to note that our pre- 
vious paper on Morse clustergE3 foreshadowed this dis- 
covery by recognizing the especial stability of the 75- 
atom Marks decahedron, thus again showing the utility 
of Morse clusters as a model system. 



C. Close-packed clusters 

The close-packed global minima are illustrated in Fig. 
|T^ . They have a diverse range of structures: there are 4 
minima that are fee, 8 that are hep and 46 that involve 
a mixture of stacking sequences and twin planes. The 
preference for close-packed structures with twin planes, 
even though at many of the sizes there are fee isomers 
with the same number of nearest neighbours, occurs for 
the same reason that bulk hep has a lower energy than 



fee for pair potentials, namely a larger energetic contri- 
bution from next-nearest neighbours. The global minima 
are broadly based on five structures which are especially 
stable: the hep 26C, the truncated octahedron 38D, the 
tetrahedral 59E and the 'twinned truncated octahedra' 
50D and 79F. The latter four give rise to peaks in A2E 
at po=U (Fig. g). 

The 38-atom truncated octahedron, 38D, is the most 
stable fee cluster in the size range we consider here. It 
becomes the global minimum at the lowest value of po 
(4.76) of any of the close-packed structures. Curiously, 
there are two ranges of po for which it is the global min- 
imum. At long range, Ennn represents a significant part 
of the total energy. The truncated octahedron is most 
stable for 4.76 < po < 5.40 because it is approximately 
spherical and so has a larger value of Ennn than the more 
oblate icosahedral structure 38E. For shorter-ranged po- 
tentials, the contribution of diminishes and so 38E 
becomes the global minimum for 5.40 < po < 6.95 be- 
cause it has a larger rt„„. Then for po > 6.95 the trun- 
cated octahedron again becomes the global minimum be- 
cause it has a lower strain energy than 38E. There is a 
growing body of experimental evidence for the impor- 
tance of truncated octahedra. Parks et al. have recently 
assigned this structure to Niag by probing the cluster's 
chemical reaetivity.EJ EXAFS (extended x-ray absorp- 
tion fine structure) spectra of small gold clusters have 
been interpreted in terms of the presence of truncated 
octahedral clusters, particularly the 38-atom truncated 
octahedronO Gold clusters passivated by alkylthiolate 
molecules selectively form truncated octahed|La|_which 
can be isolated and formed into superlattices.tilo This 
structureJs also observed for ligated 38-atom platinum 
clusters.Ea 

Structures 50D and 79F both have D^h symmetry and 
a single twin plane passing through the structure. The 
two halves of the structure have the surface morphol- 
ogy of truncated octahedra. Indeed 79F can be formed 
from the 79-atom truncated octahedron by rotation of 
one half of the structure by 60° about an axis perpendic- 
ular to one of the {111} planes. Both structures have the 
same number of nearest neighbours and 79F is slightly 
lower in energy only because of the larger contribution 
from next-nearest neighbours that results from the twin 
plane. Again there is recent experimental evidence fpr. 
the stability of this type of structure for gold clusters,a 
although at larger sizes {N=225 and 459) than consid- 
ered in this study. 

The closed-packed structures from N— 57-60 are based 
on 59E. This structure is a 31-atom truncated tetrahe- 
dron with the faces covered by four seven-atom hexagonal 
overlayers occupying hep sites with respect to the under- 
lying tetrahedron. The stability of 59E comes from the 
combination of its high proportion of {111} faces and its 
spherical shape. If atoms are added to one of the grooves 
in 59E a decahedral-like axis results, indicating a possible 
path between decahedral and closed-packed clusters. 
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FIG. 13. Global minima based upon a decahedron with five atoms along the fivefold axis. 
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FIG. 14. Global minima based upon close-packing. 
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It is worth noting that the cuboctahedron is not 
the lowest energy close-packed structure for N=55. In 
fact, it has four fewer nearest neighbours than struc- 
ture 55D. Hence, when magic numbers occur at sizes 
corresponding tn both complete Mackay icosahedra and 
cuboctahedraE^cB {N=13, 55, 147, ...) it is more likely 
that they are due to icosahedra. Furthermore, one has 
to interpret with caution those studies which seek to find 
the relative stability of fee and icosahedral structures 
by comparing cuboctahedra with Mackay icosahedraEj EJ 
because the cuboctahedra are likely to be suboptimal fee 
structures. 

Although this conclusion does not simply carry over 
to clusters surrounded by ligand shells — the ligands could 
significantly modify the relative surface energies of {111} 
and {100} faces — it is interesting to note that a recent 
reinvestigation of clusters which were origmalljy thought 
to be cuboctahedral 55-atom gold ciustersEflEj seems to 
disprove this structural assignment. E3 



D. Structures corresponding to long range 



In this subsection, we consider those structures which 
become the global minimum only at low values of po- We 
have restricted our study to those clusters with po > 3, 
since we do not know of a case where longer-ranged po- 
tentials might be relevant. The low po structures that 
we have found separate into two size ranges: those with 
A'' around 13 and those with A > 45. Xhe former set 
have been described in a prev ious paper ,£3 but we illus- 
trate them again isiHg- IIIDI because the connection to 
Kasper polyhedraEJ'Q was not originally identified and 
because they are important for understanding the struc- 
tures that occur at larger N. 

The majority of the structures associated with low po 
are polytetrahedral: the entire cluster can be divided into 
tetrahedra with atoms at the vertices. These minima 
also tend to be close to spherical in shape and highly 
strained in order to maximize n„„. As we noted ear- 
lier, the 13-atom icosahedron is polytetrahedral. In this 
case, each nearest-neighbour contact between the cen- 
tre and a vertex is the common edge of five tetrahedra. 
This is also true of all the nearest-neighbour contacts in 
the rhombic tricontahedron (45A) which do not lie on 
the surface. Nearest-neighbour contacts which are sur- 
rounded by more or less than five tetrahedra are said to 
have defects called disclination lines running along the 
interatom vector. Those contacts surrounded by more 
than five tetrahedra are termed negative disclinations (if 
there are six it is a —72° disclination, if there are seven a 
— 144° disclination, . . . ) and those surrounded by fewer 
than five tetrahedra are termed positive disclinations (if 
there are four it is a -1-72° disclination and if there are 
three a -1-144° disclination). 

Most of the structures associated with low po are poly- 
tetrahedral and involve disclinations. Although pack- 



ing five tetrahedra around a nearest-neighbour contact 
involves some strain, the energetic penalty associated 
with either more or less tetrahedra is greater. There- 
fore, structures involving disclinations are only likely to 
be lowest in energy for long-ranged potentials where the 
associated strain can be most easily accommodated, and 
where they must have a larger n„„ than the alternative 
disclination-free structures. A —72° disclination line in- 
volves less strain than a -|-72° disclination line because of 
the gap that remains when five regular tetrahedra share 
a common edge (Fig. II A (a)). Consequently, structures 
which involve only negative disclinations, or an excess 
of them, are more common amongst the low po global 
minima. 

To visualize the network of disclination lines in a struc- 
ture, one must first partition space according to the 
Voronoi procedure, in which each point is assigned to 
the Voronoi polyhedron of the atom to which it is clos- 
est. This allows nearest neighbours to be defined as those 
atoms whose Voronoi polyhedra share a face. The De- 
launay network that results from joining all such near- 
est neighbours is the dual of the Voronoi construction 
and divides all space into tetrahedra. This definition of 
a nearest neighbour has been termed geometric, rather 
than physical (e.g. using a cutoff distance), and the cor- 
responding division of space into tetrahedra is artificial 
in the sense that it is independent of whether a polyte- 
trahedral description is appropriate. 

In practice we determined the Voronoi polyhedra from 
the fact that a set of four atoms constitutes a Delau- 
nay tetrahedron if the sphere touching all four atoms 
contains no other atomsO The centre of this sphere is 
then a vertex of the Voronoi polyhedron of each atom. 
As the number of tetrahedra around a nearest-neighbour 
contact is the same as the number of sides for the face 
common to the Voronoi polyhedra of both nearest neigh- 
bours, disclination lines can then be assigned. However, 
problems can occur in assigning the Delaunay network 
if there are more than four atoms exactly on the surface 
of the sphere. Such a degeneracy, which only occurs as 
a result of symmetry, renders the analysis non-unique. 
This is the case for the 55-atom Mackay icosahedron and 
for bulk close-packed solids because of the presence of 
octahedral interstices, but it is not a problem here since 
we only apply the method to clusters that are polyte- 
trahedral in character. One further consideration is that 
the analysis should not be applied to nearest-neighbour 
contacts between the surface atoms of a cluster. 

The smallest global minimum that involves a disclina- 
tion line is 11 A where the central atom is surrounded by 
a 10-atom coordination shell. This encapsulation gives 
the structure a larger n„„ than the incomplete icosahe- 
dron (structum IIB of ref. ^8|), but results in a larger 
strain energy.tj The structure involves a single positive 
disclination line running through the centre of the cluster. 
Similarly, for A=12 and 14-16, clusters with a single co- 
ordination shell become lower in energy than structures 
based on the icosahedron at long range. Structures 11 A, 
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12A, 14A, 15A, 16A and 17F (the second lowest energy 
structure of M17 at po = 3) correspond to Kasper polyhe- 
dra. Of these structures 15 A is the most stable in terrns 
of the range of po for which it is the global minimum. E3 
The Kasper polyhedra are the deltahedral coordination 
shells that involve the minimum number of disclinations. 
They are important in the Frank-Kasper phases J^'E^ 
which are crystalline structures that are polytetrahedral 
and involve ordered arrays of negative disclination lines. 
Much interest has been focussed on the Frank-Kasper 
phases because they are closely related to icosahedral 
quasicrystals.Ej Indeed a recent three-dimensional model 
of quasicrystalline structure was based upon clusters ie-, 
volving disclination lines similar to those we find here.Ell 




FIG. 15. Low energy minima for small clusters at low val- 
ues of Po- If appropriate, the disclination network is dis- 
played next to each structure. —72° disclinations are repre- 
sented by red lines and -1-72° disclinations by blue lines. 17F 
and 20c are never global minima. 

Structures llA, 12A, 14A, 15A and 16A are all delta- 
hedral, and so growth c an oc cur in both 'anti-Mackay' 
and 'Mackay' sites (Fig. [II A ). However, if one consid- 
ers the addition of a hexagonal pyramidal cap to the 
most stable of the Kasper polyhedra, 15A, the result 
is a Dqii structure with a positive disclination running 
the length of the symmetry axis (a disclinated equiva- 
lent of the double icosahedron 19A) which has only the 
same as the icosahedral structure 22A. Disclinated 
polytetrahedra cannot compete with the disclination-free 
polytetrahedral structures produced by an anti-Mackay 
overlayer on the icosahedron. Only once the latter growth 
sequence is completed at 7V=45 are structures with discli- 
nations again global minima (Fig. 16). The one excep- 
tion is structure 38A, which is similar to the icosahedral 
structures 38B and 38C, but has two positive disclina- 
tions running through the structure in a strange double 
helical twist. 



Interestingly, structures with anti-Mackay growth on 
llA have been recently observed for M+RGat (M — metal, 
RG — rare gas) clusters where the metal ion is sufficiently 
small with respect to the rare gas atoms.Ej The metal ion 
presumably lies in the centre of the cluster, and the size 
ratio ensures that the cluster is least strained when the 
metal ion is surrounded by 10 rare gas atoms. Similarly, 
one might expect structures with negative disclinations 
to occur for AB^r clusters when A is sufficiently large 
with respect to B that a coordination number of larger 
than twelve is favoured for the A atom. 

The most stable disclinated pol ytet rahedral structures 
occur at A^=53, 57, and 61 (Fig. |l|(a)). These are the 
sizes for which complete 'anti-Mackay' overlayers on 15A, 
16A and 17F are possible. The effect of the overlayers is 
to extend the disclination lines emanating from the cen- 
tral atom. In these structures, those interior atoms not 
lying on a disclination line are icosahedrally coordinated. 
Many other minima are related to these stable structures: 
51A and 52A are based on 53A but with missing vertex 
atoms; similarly, 59B and 60A are based on 61A; and 
70A and 74A are based on 57A and 61A, respectively, 
but with an additional 13-atom cap which extends one 
of the 'arms' of the disclination network. Furthermore, 
many of the other structures include parts of the disclina- 
tion networks of 53A, 57A and 61A, but combined with 
a more disordered array of disclinations in another part 
of the cluster. 

Other interesting structures are also seen. 47A, 50A 
and 59A seem to have a mixture of Mackay and anti- 
Mackay overlayers. 64A is formed from the rhombic tri- 
contahedron 45A by the extension of the structure along 
a three-fold axis and the addition of a ring of atoms in 
the centre. In the middle of 64A is the D^d structure 
20D (the third lowest energy structure for M20 at po—5) 
which can be regarded as two highly strained interpen- 
etrating icosahedra. However, for the larger clusters it 
becomes difficult to recognize any structural motifs, and 
some just seem to be disordered tangles of disclinations. 

The above results are particularly interesting because 
of their relevance to our understanding of liquid struc- 
ture. The minima described above, when reoptimized at 
larger values of po: correlate with structures which lie in 
the lower energy range of the band of minima associated 
with liquid-like clusters. This is because siiHple liq- 
uids have significant polytetrahedral character J23 as has 
been shown by th e s u ccess of the dense random pacfc^ 
ing of hard spheresoca as a modelJor metallic glassesE3 
and later by computer simulations.EH Indeed, Nelson has 
suggested that simple liquids are polytetrahedral pack- 
ings that are charaatepzed by a disordered arrangement 
of disclination lines.OO Consequently, by examining the 
global minima associated with low po we can study the 
size evolution of polytetrahedral packings, and the devel- 
opment of bulk liquid structure. 
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At small sizes disclination-free polytetrahedra based on 
icosahedra with anti-Mackay overlayers are possible; at 
larger sizes polytetrahedra with ordered arrangements of 
disclinations are most common, and finally at the largest 
sizes in this study the polytetrahedra have a disordered 
disclination network. The latter geometries are struc- 
turally very similar to fragments of bulk liouid, except 
that the density of disclination lines is lower .EJ 

Theoretical studies of sodium clusters have shown that 
amorphous structures are lower in energy than regular ge- 
ometries up to at least 340 atoms, the largest size consid- 
ered in that study.ll23 The present results suggest that the 
these disordered structures are due to the relatively long 
range of the sodium interatomic potential. Amorphous 
structure for sodium is also suggested by the experim 
tal observation of electronic shellsE21l and supershells 
which are incompatible with— any of the ordered mor- 
phologies one might expect.E£3 Indeed, the transition 
from electronic to geometcic magic numbers which oc- 
curs at about 1000 atomsc3 probably reflects a change in 
the lowest energy structures from amorphous minima to 
Mackay icosahedra, which can be explained by tiif|-t,win 
effects of the size and the range of the potentialOH 

Since most minima associated with small values of po 
do not have a common lattice or packing scheme, it is dif- 
ficult to predict the low energy structures, and so most 
have been found by one of the global optimization meth- 
ods. Fortunately, these methods are most likely to suc- 
ceed for long range potentials because the correspond- 
ing PES's are smoother and support fewer minima than 
for shorter range. Another consequence of the lack of 
lattice structure is that the differentiation between near- 
est neighbours and next-nearest neighbours becomes am- 
biguous and the decomposition of cqn. (2) more arbitrary. 
For these clusters we chose to define n„„ using a cutoff 
at po—4: which corresponds to the distance at which the 
pair energy is 0.6 e. 



IV. DISCUSSION 

In this paper we have attempted to find the global min- 
ima for Morse clusters as a function of po and the number 
of atoms. The global potential energy minimum repre- 
sents the equilibrium structure at zero Kelvin, but to 
predict the structure at non-zero temperatures we must 
consider the free energy, and the effect of other low en- 
ergy minima. This can be done by summing the density 
of states over all the relevant minimaJi£j An illustration 
of this approach for understanding tbfi.|Sttucture of Morse 
clusters has been given previously,HE23 revealing that 
for M75 at po=6, the equilibrium structure changed from 
decahedral to icosahedral at very low temperature. This 
transition is simply a consequence of the larger entropy 
for the icosahedral region of configuration space — there 
are far more low energy icosahedral minima. 

This effect is likely to be general. At the magic num- 



bers for a morphology the entropy will be low because 
there is a single unique low entropy structure and a 
large gap to other isomers with the same morphology. 
Therefore the finite temperature equivalent of the struc- 
tural phase diagram of Fig. |^ is likely to show weakened 
magic number effects and so have smoother boundaries 
between the different morphologies. Furthermore, the 
energy gap between the lowest energy ordered minimum 
and the liquid-like band of minima increases as the range 
becomes shorter_-and hence the melting temperature in- 
creases with po-t3 Therefore the region of the phase di- 
agram where disordered polytetrahedral structures have 
the lowest free energy is likely to spread up from the 
bottom as the temperature increases. 

In this paper we have considered only isotropic pair- 
wise additive interactions. As noble gas clusters and clus- 
ters of Ceo molecules can be reasonably modelled by such 
potentials we would expect the structures we have found 
at the appropriate values of po to be similar to the ac- 
tual structures of these clusters. Our results lead us to 
predict that neutral clusters of Ceo molecules exhibit dec- 
ahedral and fee structures at small sizes because of the 
short range of the intermolecular potential. This basic 
conclusion has beeit-eonfirmed in studies using more re- 
alistic potentials OE2a 

In contrast, making predictions for metal clusters is 
problematic because the range of the potential is only one 
factor determining the structure and, many - body terms. 
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These terms 



in particular, may also be important, 
may affect the relative surface energies of {111} and 
{100} faces, and so alter the energetic competitiaa be- 
tween icosahedral, decahedral and fee structures.Eil For 
example, in a study of lead clusters cuboctahedra are al- 
ways found to be lower in energy than icosahedra because 
the susfece energies of {111} and {100} faces are nearly 
equal .H 

Nevertheless our results are of value to the field of 
metal cluster structure. Firstly, they enable particularly 
stable structural forms to be identified. For example, in 
our previous paper on Morse clusters we identified the 38- 
atom octahedron and the 75-atom Marks decahedron as 
particularly stable. Subseojoently, they have both been 
observed experimentally;EilLJ it even being possible to 
isolate fractions of the latter when passivated by surfac- 
tants. This correspondence between the Morse structures 
and those of real systems encourages us to believe that 
some of the general principles that determine stability in 
our simple model system do carry over to real clusters. 

Secondly, the Morse structural database should be use- 
ful in providing candidate structures for comparison with 
the indirect structural information yielded by experi- 
ments on size-selected clusters. Finally, the database can 
also provide plausible starting structures for theoretical 
studies with more realistic, but computationally expen- 
sive, descriptions of the interactions; this expense would 
prevent the type of extensive searches that have been per- 
formed in this paper. Indeed, we have used the database 
in this way in studies of metal clusters modelled by the 
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Sutton-Chejuiamily of potentialscJ and clusters of Cg 
molecules JI3'E23 For these reasons, the coordinates fc 
the global minima given in this and previous naaerf'^^''*^' 
will be made available on the world-wide-web 
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V. CONCLUSION 

We have shown how the range of an isotropic pairwise 
additive potential determines the structure of atomic 
clusters. In particular, we have identified four principal 
structural regimes. For long-ranged potentials at sizes 
TV = 10-18 and > 45 the global minima are generally 
polytetrahedral and can be analyzed in terms of discli- 
nation lines. For the smaller clusters the disclinations 
form ordered arrangements and the structures are frag- 
ments of bulk Frank-Kasper phases. As the size increases 
it becomes more likely that the disclination network is a 



disordered tangle; these global minima have amorphous 
structures similar to those exhibited by liquid-like clus- 
ters. At intermediate ranges of the potential icosahedra 
are dominant. As the range decreases, first decahedral 
and then fee structures become lowest in energy. These 
trends have been explained by considering the strain en- 
ergies and the number of nearest neighbour contacts as- 
sociated with each regime. The effect of decreasing the 
range of the potential is to destabilize the strained struc- 
tures. 
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TABLE I. Global minima of M]v for TV < 80. Energies at values of po for which the structure is the global minimum are 
given in bold, pmin and pmax give the range of po for which the structure is the global minimum. If at a particular value of 
Po a structure is not a minimum but a higher order saddle point, the index of the stationary point (the number of negative 
eigenvalues of the Hessian) is given in square brackets after the energy. -Bstrain has been calculated at po=10. If a structure is 
not stable at po = 10 no value of -Bstrain is given. All energies are given in e. 





PG 




-^'strain 


Po = 3 


po = 6 


po = 10 


po = 14 


pmin 


pmax 


5A 


Dsn 


9 


0.000 


-9.299500 


-9.044930 


-9.003565 


-9.000283 






6A 


Oh 


12 


0.000 


-13.544229 


-12.487810 


-12.094943 


-12.018170 






7A 


D,,h 


16 


0.062 


-17.552961 


-16.207580 


-15.956512 


-15.883113 






8A 


Did 


18 


0.006 


-22.042901 


-19.161862 


-18.275118 


-18.076248 




5.28 


SB 


Cs 


19 


0.062 




-19.327420 


-18.964638 


-18.883688 


5.28 




9A 


Dsn 


21 


0.002 


-26.778449 


-22.330837 


-21.213531 


-21.037957 




3.42 


9B 


Civ 


23 


0.186 


-26.607698 


-23.417190 


-22.850758 


-22.644892 


3.42 




lOA 


Did 


24 


0.002 


-31.519768 


-25.503904 


-24.204958 


-24.031994 




2.28 


lOB 


Czv 


27 


0.694 


-31.888630 


-27.473283 


-26.583857 


-26.132735 


2.28 




llA 


Did 


34 


10.374 


-37.930817 


-28.795153(4] 


-23.666072(5] 






3.40 


UB 


Civ 


32 




-37.891674[1] 








3.40 


3.67 


lie 


Civ 


31 


0.792 




-31.521880 


-30.265230 


-29.588130(1] 


3.67 


13.57 


UD 


Cs 


30 










-29.596054 


13.57 


15.29 


HE 


Ci 


30 


0.248 


-36.613190(1] 


-30.698890 


-29.808994 


-29.524398 


15.29 


20.60 


UF 


Civ 


29 


0.001 


-36.697760 


-30.431713 


-29.215510 


-29.037941 


20.60 




12A 


Civ 


38 




-43.971339(1] 


-35.199881(1] 








2.63 


12B 


Csv 


36 


1.704 


-44.097880 


-36.400278 


-34.366755 


-33.115942(1] 


2.63 


12.15 


12C 


Cs 


34 










-33.199505 


12.15 


13.03 


12D 


Did 


34 


0.346 


-41.816393 


-34.838761 


-33.724155 


-33.332305 


13.03 


17.08 


12E 


Dzh 


33 


0.001 


-42.121440 


-34.568002 


-33.222331 


-33.038298 


17.08 




13A 


h 


42 


2.425 


-51.737046 


-42.439863 


-39.662975 


-37.258877 




14.76 


13B 


Dsh 


37 


0.141 


-49.998058(1] 


-39.360710(1] 


-37.208019(1] 


-36.790507 


14.76 




14A 


Civ 


46 


4.373 


-56.754744 


-44.827522(1] 


-41.717043(1] 






3.23 


14B 


Czv 


45 


2.425 


-56.660471 


-45.619277 


-42.675222 


-40.259823 


3.23 


13.06 


14C 


Civ 


41 


0.141 


-55.971620(2] 


-43.634048(1] 


-41.249282 


-40.798348 


13.06 




15A 


Ded 


50 


9.527 


-63.162119 


-47.570579 


-40.569211(10] 


-35.758904(13] 




3.72 


15B 


Civ 


49 


2.573 


-62.593904 


-49.748409 


-46.541404 


-44.086633 


3.72 


12.71 


150 


Civ 


45 


0.141 


-62.631372(2] 


-47.952559 


-45.293844 


-44.806437 


12.71 




16A 


Dzh 


54 


11.222 


-69.140648 


-50.834213 


-42.887569(12] 






3.39 


16B 


Cs 


53 


2.868 


-68.757203 


-53.845835 


-50.261947 


-47.831957 


3.39 


11.99 


16C 


Civ 


49 


0.142 


-68.575718(1] 


-52.265348 


-49.338173 


-48.814517 


11.99 





18 



TABLE I. continued. 





PG 


n 


'Strain 


po = 3 


po = 6 


Po = 10 


Po = 14 


Pmin 


f-'TTlClX 


17A 


^ on. 


58 




-75.662417 


-53. 156042 [2] 










3.42 


17B 


Czv 


57 


3.372 


-75 147372 


-57 884517 


-53 772213 


-51 329560 


3 


,42 


4.88 


17C 


Cs 


57 


3.281 


-75.091367 


-57.912963 


-53.862044 


-51.440588 


4 


.88 


4.91 


17D 


Co 


57 


3.163 


-75 005403 


-57.941386 


-53 983559 




4 


.91 


11.30 


17E 




53 


0.142 


-74 868921 fll 


-56 573571 

XJ\J • 'J 1 'J 'J 1 i- 


-.53 382277 


-52 822.'>88 

\J ■ Q__J 


11 


.30 




18 A 


C2 


65 




-82.579266 


-59.881449[1] 






2, 


.14 


3.03 


18B 




62 


4.500 


-82 548885 


-62.689245 


-57 657135 


-54 059707 


3 


.03 


10.13 


18C 


^ s 


61 


3.528 


-81 256639 


-62.002920 


-57 634324 


-55 776126 


10 


.13 


10.43 


18D 


Drh 


57 


0.142 


-81.490185 


-60.926500 


-57.429683 


-56.830907 


10, 


.43 




19A 


Drh 


68 


6.001 


-90.647461 


-68.492285 


-62.166843 


-56 676685f4l 






10.70 


19B 




61 


0.151 


-87.485744 


-65.064771 


-61.427105 


-60.812425 


10 


.70 




20A 


Ch 


72 


6.507 


-97.417393 


-72.507782 


-65.679115 


-61.327229 


2, 


.02 


10.24 


20B 


Co,. 


65 


0.162 


-94.222416 


-69 202704 


-65 423697 


-64.791953 


10, 


.24 




21A 




76 


7.133 


-104.336946 


-76.487266 


-69.068687 


-65.179591[1] 






5.40 


21B 




76 


6.922 


-104.004129 


-76 529139 


-69 276346 


-65 778898 


5, 


.40 


9.84 


21C 




69 


0.169 




-73 577014 


-69 449904 


-68 783571 


9, 


.84 




22A 


Cs 


81 


8.198 


-112.041223 


-81.136735 


-73.014321 


-68.580862[1] 






9.58 


22B 


C^ 


73 


0.169 




-77 887855 


-73.494292 


-72.791747 


9, 


.58 




23A 




87 


10.597 


-120 786879 


-86 735494 


-76 630624 


-70 816059f5l 






8.35 


23B 




84 


6.097 


-116.279438 


-84.940552 


-78.143867 


-74.442341 [2] 


8, 


.35 


9.71 


23C 


c. 


78 


0.433 




-83.504908 


-78.325380 


-77.302495 


9, 


.71 


22.45 


23D 




76 


0.007 




-82.252747 


-76 898804 


-76 157457 


22 


.45 




24A 


Cs 


91 


10.947 


-127.884549 


-90.685398 


-80.295459[1] 


-74.613738(2] 


2 


.79 


8.55 


24B 


Ci 


82 


0.434 




-87.820376 


-82.370214 


-81.309508 


8, 


.55 


15.34 


24C 




81 


0.007 




-87 626843 


-81.941420 


-81 164168 


15, 


.34 




25A 


Cs 


96 


12.090 


-136.072704 


-95.127899 


-84.168765 




2, 


.62 


7.76 


25B 




87 


0.893 




-93.342771 


-86.989688 


-85.477376 


7, 


.76 


15.40 


25C 


c. 

^ s 


85 


0.010 




-92.241466 


-86 015371 


-85 176789 


15, 


.40 




26A 


a 


102 


15.393 


-145.322134 


-100.549598 


-86.882333[11] 








7.89 


26B 




91 


0.491 


-138.940920[3] 


-97.363225 


-91.370250 


-90.210764 


7, 


.89 


14.20 


26C 


J-^ oil 


90 


0.008 




-97 648652 


-91 085136 


-90 189274 


14, 


.20 




27 A 


Cs 


106 




-152.513867 


-104.489430 










4.72 


27B 




106 




-151.734925 


-104.745275 






4, 


.72 


7.71 


27C 


c. 


97 


3.076 




-102.749592 


-94 736230 


-91 236563fll 


7, 


.71 


8.53 


27D 


Cs 


95 


0.491 


-146.171281 [2] 


-101.722918 


-95.419490 


-94.219798 


8, 


.53 


14.22 


27E 


Cs 


94 


0.008 




-101.920210 


-95.121297 


-94.195553 


14, 


.22 




28A 


c. 

^ s 


111 




-160 773356 


-108 854564 










4.82 


28B 


Cs 


111 




-160.385239 


-108.997831 






4, 


.82 


6.70 


28C 




102 


3.369 




-108.186446 


-99.524026 


-95.692056 


6, 


.70 


8.88 


28D 




100 


0.983 




-107 213896 


-100 008358 


-98 331711 


8, 


.88 


14.57 


28E 




98 


0.008 




-106.238844 


-99.161073 


-98.202117 


14, 


.57 




29A 




117 




-170 115560 


-114.145949 










6.90 


29B 


Ci 


106 


3.370 




-112.655980 


-103.586316 


-99.702681 


6, 


.90 


8.95 


29C 


Ci 


104 


0.988 




-111.543685 


-104.051412 


-102.333194 


8, 


.95 


11.07 


29D 




103 


0.246 




-111.508961 


-103.946446 


-102.743899 


1 1 


.u / 


11.17 


29E 




103 


0.225 


-164.720567 


-111.353973 


-103.926827 


-102.774589 


11, 


.17 


21.90 


29F 


Czv 


102 


0.010 




-111.135014[1] 


-103.298583 


-102.227015 


21, 


.90 




30A 


C2V 


121 




-177.578647 


-118.115802 










4.47 


30B 


C2v 


121 


8.279 


-176.940971 


-118.432844 


-103.407917[1] 




4 


.47 


6.86 


30C 


C2v 


109 


1.538 




-117.010672 


-108.571179 


-106.426033 


6, 


.86 


12.61 


30D 


Cs 


107 


0.238 




-115.940968 


-107.985497 


-106.765372 


12, 


.61 


21.33 


30E 


Cs 


106 


0.012 




-115.625207 


-107.367696 


-106.239394 


21, 


.33 
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TABLE I. continued. 





PG 


n 


-^/strain 


Po ~ 3 


Po — 6 


Po = 10 


Po — 14 


Pmin 


771 a. X 


31A 


Ci 


126 






-185.984248 


-122.440052 










4.60 


31B 


C. 


126 






-185.299446 


-122 857743 






4, 


.60 


6.35 


31C 




115 


3 


.490 




-122.342421 


-112.542723 


-108.722594 


6, 


,35 


8.93 


31D 




112 


0, 


.269 


-184.500261 


-121.523693 


-113.066115 


-111.724317 


8, 


,93 


10.98 


31E 




112 


0, 


.246 


-181.901112[1] 


-121.367441 


-113 048378 


-111.760670 


10, 


,98 


20.90 


31F 


c. 


111 


0, 


.011 




-120.805231 


-112.376045 


-111.239641 


20, 


,90 




32A 


C'2v 


132 






-195 468461 


-127 643751 










5.93 


32B 




120 


3 


.831 




-127.771395 


-117.284334 




5, 


,93 


10.27 


32C 


Cs 


116 


0, 


.269 




-125.953587 


-117.116756 


-115.732921 


10, 


,27 


11.33 


32D 




116 


0, 


.260 




-125 950348 


-117 105640 


-115 748493 


11, 


,33 


11.40 


32E 


^ s 


116 





.246 


-188.807956[1] 


-125.644966 


-117.086755 


-115.767561 


11, 


,40 


20.90 


32F 


Cs 


115 


0, 


.011 




-125.126341 


-116.415892 


-115.246208 


20, 


,90 




33A 


^ s 


137 






-204.208737 


-131 704206fll 










5.35 


33B 




137 






-203.575130 


-131.773513 






5, 


,35 


5.74 


33C 


^ s 


124 


3 


.832 




-132.287431 


-121.352047 




5, 


.74 


8.46 


33D 




121 


0, 


.313 




-131 555811 


-122.160387 


-120 665738 


8, 


.46 


9.71 


33E 




121 


0, 


.269 




-131.378662 


-122.167479 


-120.741345 


9, 


.71 


19.88 


33F 


^ s 


120 


0, 


.011 




-130 466899 


-121.454645 


-120.252542 


19. 


,88 




34A 




143 






-214.068392 


-136 468311 










5.87 


34B 




128 


3, 


.832 


-209.946801 


-136.797544 


-125.419532 




5, 


.87 


8.52 


34C 




125 


0, 


.313 




-135.989024 


-126.211991 


-124.674640 


8, 


.52 


9.55 


34D 




125 





.289 




-135 963257 


-126.217685 


-124 717107 


9, 


.55 


10.72 


34E 


Cs 


125 


0, 


.269 


-206.153688(2] 


-135.656902 


-126.205973 


-124.748271 


10. 


,72 


19.86 


34F 




124 


0, 


.013 




-135.532334 


-125.625477 


-124.283737 


19. 


,86 




35A 


c. 


147 






-221.771452 


-140 503355 










4.05 


35B 




147 






-221.293580 


-141.106305 






4, 


,05 


5.62 


35C 


Ci 


133 


3 


.975 




-141.957188 


-130.282153 


-126.119817 


5, 


,62 


8.08 


35D 




130 


0, 


.304 




-141.402997 


-131.274427 


-129.699893 


8, 


,08 


11.46 


35E 




130 


0, 


.279 


-215.495670[2] 


-141.051852 


-131.239157 


-129.737360 


11, 


,46 


19.59 


35F 




129 


0, 


.011 




-140.132124 


-130.533259 


-129.265442 


19, 


,59 




36A 


^ s 


152 






-230.508264 


-144.827464 










4.68 


36B 


Ci 


152 






-229.689435 


-145.273702 






4, 


,68 


5.07 


36C 


Cs 


138 


4 


.351 




-147.381965 


-134.989116 




5, 


,07 


9.51 


36D 


Ci 


134 


0, 


.334 




-145 984498 


-135 312738 


-133 656022 


9, 


.51 


9.86 


36E 


c. 


134 





.303 




-145.732293 


-135.315392 


-133.706961 


9, 


.86 


11.41 


36F 


^ s 


134 


0, 


.280 


-222 900558fll 


-145 378315 


-135.281551 


-133.744666 


11, 


,41 


19.56 


36G 


Cs 


133 


0, 


.015 




-145.649065 


-134.766521 


-133.308212 


19, 


,56 




37 A 


^ s 


158 






-240 008130 


-149 765993 










5.25 


37B 


Ci 


142 


4, 


.352 




-151.891203 


-139 055858 


-134 764790 


5, 


.25 


7.62 


37C 


C2v 


139 





.357 




-151.435570 


-140.360228 


-138.622227 


7, 


.62 


10.17 


37D 


C'Zv 


139 





.306 


-233.080175(3] 


-151.061897 


-140.354330 


-138.708582 


10, 


,17 


18.49 


37E 




138 


0, 


.015 




-151.100849 


-139.837570 


-138.320817 


18, 


,49 




38A 


Do 


164 


30, 


.138 


-249.159174 


-153 208710f3l 


-134 319674f3l 








3.15 


38B 


Ci 


163 


27, 


.962 


-249.087740 


-154.165069 


-135.519468 


-129.339213 


3, 


.15 


4.70 


38C 


c. 


163 






-248.600369 


-154.041575 






4, 


.70 


4.76 


38D 


Oh 


144 


0, 


.013 


-246.414723(4] 


-157 406902 


-145 849817 


-144.321054 


4, 


,76 


5.40 


38E 




147 


3 


.714 




-157.477108 


-144.756555 




6, 
5, 


.95 
.40 


6.95 


39A 


C'Zv 


169 






-258.944962 


-158 266473 










4.48 


39B 


Cbv 


153 


5, 


.039 




-163.481990 


-149.455758 




4, 


,48 


9.46 


39C 


C'iv 


148 


0, 


.014 




-161.701728 


-149.887173 


-148.327400 


9, 


.46 




40A 


Cs 


174 






-268.394773 


-163.956596 










4.51 


40B 


Cs 


157 


5, 


.040 


-263.925318 


-167.993097 


-153.524033 


-148.492778(1] 


4, 


.51 


9.51 


40C 


D4h 


152 


0, 


.014 


-261.115599(3] 


-165.996196 


-153.924517 


-152.333745 


9, 


.51 




41A 


Czv 


180 






-278.405573 












4.60 


41B 


('.. 






,()1() 


-27:!.()99r,()«[i] 


-172.526828 


-i.")7.r,<):!2()() 




1 


,()() 


.S.75 
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TABLE I. continued 





PG 


n 


-^strain 


Po — 3 


po ~ 6 


po = 10 


po = 14 


Pmin 


Pmax 


41C 




157 


0, 


.435 




-171.124667 


-158.505557 


-156.511959 


8. 


.75 


9 


.60 


41D 




157 


0, 


.368 


-269 034859 f41 


-170 808262 


-158 520688 


-156 630647 


9, 


.60 


11, 


.06 


41E 




157 


0, 


,366 


-267.688391(4] 


-170.758682 


-158.518482 


-156.633479 


11, 


.06 


16, 


.83 


41F 


^ s 


156 


0, 


.017 




-170 902496 


-158 066916 

A.\J\^ *\J\J\J \J \.\J 


-156 359354 


16, 


.83 






42A 




186 






-288.335415 


-172 252388fll 










4 


.54 


42B 


c, 

^ s 


166 


5, 


.275 


-282.296210[1] 


-177.680222 


-162.365546 


-157.124535 


4, 


.54 


9 


.77 


42C 


Cs 


161 


0, 


.368 




-175.184389 


-162.565848 


-160.638215 


9, 


.77 


11, 


.10 


42D 


^ s 


161 


0, 


.367 




-175 135897 


-162 563540 


-160.641020 


11, 


.10 


16, 


.82 


42E 


Cs 


160 


0, 


.017 




-175.405455 


-162.133633 


-160.371244 


16, 


.82 






43A 


C'2v 


192 






-298.172449 


-175 540353 f2l 










4, 


.40 


43B 


c. 


171 






-290.632163[2] 


-183.092699 






4, 


.40 


9 


.34 


43C 




166 


0, 


.443 




-180.837781 


-167.585928 


-165.511347 


9, 


.34 


9 


.57 


43D 


^ 2.V 


166 


0, 


.374 


-288 301568[2l 


-180 515098 


-167.602632 


-165.634973 


9, 


.57 


16, 


.66 


43E 


Ci 


165 


0, 


.017 




-180.662511 


-167.152124 


-165.372773 


16, 


.66 






44A 


^ ov 


198 






-308 277011 


-179 032310f8l 










4 


.38 


44B 




175 


5. 


.722 


-302 003640 


-187.626292 


-171 071131 


-165 689737 

A~\J^ »\JKJ I (J I 


4, 


.38 


9 


.38 


44C 


Cs 


170 


0, 


.443 




-185.214201 


-171.631403 


-169.519034 


9, 


.38 


9 


.65 


44D 


Cs 


170 


0, 


.374 


-300.113694(3] 


-184.841212 


-171.645465 


-169.642441 


9, 


,65 


16, 


.65 


44E 


Ci 


169 


0, 


.020 




-185.296321 


-171.226450 


-169 385406 


16, 


.65 






45A 




204 






-318 660653 


-182 301077f24l 










4, 


.24 


45B 


Ci 


180 


5. 


.965 


-310.750385 


-192.954739 


-175.389424 


-169.712668 


4, 


.24 


8, 


.72 


45C 




175 


0, 


.691 




-191 230983 


-176 548074 


-174 095718 


8, 


.72 


8 


.95 


45D 




175 


0, 


.452 




-190 548331 


-176.665997 


-174 510080 


8, 


.95 


11 


.79 


45E 




175 


0, 


.451 


-303.604189(3] 


-190.500098 


-176.663202 


-174.511633 


11, 


.79 


15, 


.03 


45F 




174 


0, 


.018 




-190 708793 


-176 296327 


-174 397893 


15, 


.03 






46A 


Cs 


207 






-327.033118 












3 


.96 


46B 


C-Zv 


186 


6, 


.646 


-320.118738(1] 


-199.177751 


-181.236182 


-174.605103(1] 


3, 


.96 


10 


.46 


46 C 


Cs 


179 


0, 


.452 




-194 923850 


-180 711434 


-178 517769 


10, 


.46 


11, 


.79 


46D 




179 


0, 


.451 




-194 876995 


-180 708654 


-178 519320 


11, 


.79 


14 


.97 


46E 


Cs 


178 


0, 


.019 




-195.426416 


-180.401063 


-178.416567 


14, 


.97 






47A 


Ci 


210 


40, 


,212 


-336.666189 


-194.242255 


-170 368382 








4 


.03 


47B 


Ci 


190 


6, 


,648 


-331.591008 


-203.704178 


-185.301500 


-178.621197 


4, 


.03 


9 


.60 


47C 


C'2v 


184 


0, 


,461 


-323 010894f6l 


-200 256505 


-185.745782 


-183 508227 


9, 


,60 


14, 


.86 


47D 


c. 


183 


0, 


.019 




-200.473614 


-185.381591 


-183.411312 


14, 


.86 






48A 




213 






-346 662788 












4, 


.01 


48B 


^ s 


195 


7, 


.464 


-340 395646fll 


-209.044000 


-189 5566723 

A-\Jj \J ^XJXJ\J\J J 


-182 617917fll 

X . X 1 ty X 1 X 


4, 


,01 


8 


.01 


48 C 


C'2v 


190 


0, 


.808 




-207.541529 


-191 586299 


-188 888965 


8, 


,01 


20, 


.22 


48D 


Cs 


187 


0. 


.020 




-205.418604 


-189.535478 


-187.440413 


20, 


.22 






49A 


c. 


219 






-356.412817 












3 


.91 


49B 




201 


7, 


.770 


-350 540378fll 


-215.253702 


-195 275059 


-187.212009 


3, 


.91 


9 


.70 


49C 


^ s 


194 


0, 


.808 




-211 978083 


-195.639320 


-192.898412 


9, 


.70 


16, 


.14 


49D 




192 


0, 


.021 




-211.117300 


-194.627748 


-192.455390 


16, 


.14 






50A 




224 






-366.635589 












4, 


.04 


SOB 


^ s 


205 


7. 


.770 




-219.820229 


-199.344795 


-191.239955 


4, 


.04 


8, 


.85 


50C 


c. 


199 


0. 


.809 


-360.600022(2] 


-217.406110 


-200.691615 


-197.906988 


8, 


.85 


10 


.50 


SOD 


J-^ on. 


198 


0. 


.020 




-217.432647 


-200 640183 


-198 455633 


10, 


.50 






51A 


C2 


232 






-376.673413 


-217.627962 










3 


.74 


51B 




210 


9, 


.360 


-373.932555(1] 


-225.391240 


-202.911593 




3, 


.74 


6 


.59 


51C 


c. 

^ s 


210 


8, 


.551 


-374 031668 1 


-225 236760 


-203 643178 


-195 286917 


6, 


.59 


9 


.18 


51D 


Cs 


203 


0, 


.809 




-221.842387 


-204.744598 


-201.916433 


9, 


.18 


10 


.33 


51E 


Cs 


202 


0. 


.023 




-222.080140 


-204.714809 


-202.468274 


10, 


.33 






52A 


Cs 


238 






-387.587332 












3 


.71 


52B 


Csv 


216 


9, 


.618 


-384.520749 


-231.615013 


-208.673261 


-199.353198 


3, 


.71 


6 


.61 


52C 


C2v 


216 


8, 


.840 




-231.443769 


-209.379230 


-199.715952 


6, 


.61 


9 


.72 


52D 


Cz. 


208 


0, 


,809 


-.380.031663[.3] 


-227.268275 


-209.796728 


-206.924957 


9, 


,72 


10 


.15 


52E 


Clr 


207 


0, 


,022 




-227.500881 


-201). 7« 1«:!2 


-207.480764 


10, 


,1.") 
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TABLE I. continued 





PG 


n 


-^strain 


po ~ 3 


Po — 6 


Po = 10 


Po = 14 


Pmin 


771 a. X 


53A 


^ba 


244 






-398.783184 












3.70 


53B 




222 


9 


.902 


-395 150063 


-237 834976 


-214.409721 


-203 180586 fll 


3 


.70 


10.30 


53C 


Cs 


211 





.025 




-232.147467 


-213.859451 


-211.493405 


10, 


.30 




54A 


^ s 


248 






-407 966010 


-232 987285 










3.41 


54B 




228 


10 


.215 


-405 831836 


-244.058174 


-220.118611 


-208 329991 


3 


.41 


10.32 


54C 


C2v 


217 





.465 




-238.114592 


-219.426691 


-216.636864 


10 


.32 


15.10 


54D 




216 





.024 




-237 568575 


-218.929453 


-216 505895 


15, 


.10 




55A 


^ i 


252 






-417.918562 












3.25 


55B 


Ih 


234 


10 


.543 


-416.625645 


-250.286609 


-225.814286 


-213.523774 


3 


.25 


11.15 


55C 




221 





.465 




-242.622450 


-223.482018 


-220.646208 


11, 


.15 


15.09 


55D 


^ s 


220 





.027 




-242.211381 


-223.003991 


-220.518533 


15, 


.09 




56A 


Cs 


258 






-428.611289 












3.60 


56B 




237 


10 


.545 


-425 709433 


-253.922955 


-228.900154 


-216 537609 


3 


.60 


10.18 


56C 


C'2v 


226 





.481 




-248 167730 


-228.552734 


-225 633470 


10, 


.18 


11.38 


56D 


Cs 


226 





.465 




-248.051289 


-228.533805 


-225.655136 


11, 


.38 


15.07 


56E 


^ oil 


225 





.026 




-247.635625 


-228.074040 


-225.531024 


15, 


.07 




57A 


^ on. 


264 






-439.960320 












4.13 


57B 


Cs 


241 


10 


.563 


-436.124573[1] 


-258.041717 


-232.877539 


-220.502856 


4 


.13 


9.61 


57C 




231 





.465 


-426 761701 3 


-253.482114 


-233 585699 


-230 663986 


9, 


.61 


21.83 


57D 


c. 


229 





.028 




-252.749395 


-232.226794 


-229.557450 


21, 


.83 




58A 




267 






-449.432282 


-253.758845(1] 










4.05 


58B 


Csv 


246 


10 


.585 


-442.501104(1] 


-263.410755 


-237.898423 


-225.469489 


4, 


.05 


9.53 


58C 


^ on 


237 


1 


.467 




-260.086805 


-238.658866 


-234.809078 


9, 


.53 


14.63 


58D 


Ci 


235 





.465 




-257.876211 


-237.631296 


-234.671473 


14, 


.63 


14.94 


58E 




234 





.027 




-258.241372 


-237 298283 


-234 569986 


14, 


.94 




59A 




267 






-459.509280 


-260.795085(1] 










3.77 


59B 




272 






-459.250400 


-261.475383(2] 






3, 


.77 


4.26 


59C 


Ci 


250 


10 


.832 




-267 857802 


-241 705877 


-229.053528 


4, 
7, 


.26 
.93 


4.86 
9.09 


59D 


C2 


250 


10 


.939 


-453.109891[2] 


-267.945226 


-241.627254 


-228.849111 


4, 


.86 


7.93 


59E 


TrI 


240 





.027 




-264.710286 


-243 330508 


-240.572493 


9, 


.09 




60A 


Cs 


278 






-470.448485 












4.51 


60B 


^ s 


255 


11 


.042 




-273.341243 


-246.579281 


-233.702420 


4, 


.51 


9.29 


60C 


C'2,1! 


246 


1, 


.470 




-269 958055 


-247 763712 


-243 823553 


9, 


.29 


9.85 


60D 


C, 


245 





.540 




-269.408035 


-247.785988 


-244.566234 


9 


.85 


11.00 


60E 




245 





.536 




-269.292014 


-247.781611 


-244.572017 


11, 


.00 


13.95 


60F 


c, 

^ s 


244 





.027 




-269 050233 


-247 370671 


-244.579066 


13, 


.95 




61A 


Trl 


284 






-482 025765 


-271.692826 










4.65 


61B 


C2v 


260 


11 


.174 




-278.726626 


-251.502374 


-238.498558 


4, 


.65 


9.54 


61C 




249 





.027 




-274.114496 


-252.402685 


-249 587740 


9, 


.54 




62A 


^ s 


293 






-491.052378 


-272.267523 










4.00 


62B 


C, 


264 


11 


.707 


-485.425587[3] 


-283.183002 


-255.038209 


-241.549494 


4, 


.00 


6.33 


62C 




264 


11 


.174 




-283.102592 


-255 558891 


-242 509563 


6, 


.33 


9.22 


62D 




255 


1 


.473 




-279.826688 


-256.868159 


-252.837765 


9, 


.22 


9.64 


62E 


Ci 


254 





.623 




-279.851878 


-256.899841 


-253.442114 


9 


.64 


9.78 


62F 


^ s 


254 





.560 




-279.460422 


-256.908198 


-253 554702 


9, 


.78 


10.47 


62G 


Ci 


254 





.526 




-279.357773 


-256.897318 


-253.612942 


10, 


.47 


14.14 


62H 


Ci 


253 





.027 




-278.499195 


-256.445284 


-253.594448 


14, 


.14 




63A 


^ s 


301 






-501.731893 












3.83 


63B 


Ci 


269 


11 


.880 


-497.181562 


-288.560948 


-259.921304 


-246.279386 


3, 


.83 


8.87 


63C 


Cs 


259 





.533 


495.356854(1] 


-285.141958 


-262.030007 


-258.620607 


8, 


.87 


14.14 


63D 


Cs 


258 





.027 




-283.876854 


-261.487847 


-258.601156 


14, 


.14 




64A 


Dsd 


307 






-512.831683 












3.85 


64B 


Cs 


274 


12 


.066 


-503.234962 


-293.931716 


-264.790492 


-250.989419 


3 


.85 


8.22 


64C 


C2V 


265 





.552 


-501.818445[4] 


-291.412250 


-268.023828 


-264.587042 


8, 


.22 


24.31 


64i) 


r. 


2()2 









-2«<).:i:i:!r)87 


-2()5.722():!o 


-2()2.()iir)i)i) 


2i 
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TABLE I. continued 





PG 


n 


-^strain 


= 3 


Po — 6 


Po = 10 


Po = 14 


pmin 


f 'Tn a. X 


65A 




310 






-523.446427 












3.86 


65B 




278 


12. 


.494 




-298 321863 


-268 406866 


-254 407364 


3 


.86 


4.74 


65C 


C2 


278 


12, 


.574 


-515.414573 


-298.392345 


-268.351634 


-254.301092 


4. 


.74 


6.32 


65D 


^ 1 


278 


12, 


.081 




-298 317501 


-268 836726 


-254.984242 


6. 


.32 


8.12 


65E 


^■^ zv 


269 


0, 


.603 




-296 31 7862 


-1T2 098786 


-268 51 3581 


8. 


.12 


10.53 


65F 


Ci 


269 


0, 


.565 




-296.027614 


-272.083782 


-268.578320 


10. 


.53 


11.30 


65G 




269 


0. 


.552 




-295 806729 


-272 069820 


-268 594702 


11. 


.30 


19.71 


65H 


^ s 


267 


0. 


.027 




-293 637761 


-270 572952 


-267 614569 


19. 


.71 




66A 


^ s 


314 






-534.464040 












3.82 


66B 




284 


15. 


.128 


-530 516639 


-303 512725 


-271.712472 


-256.391481 


3. 


.82 


4.83 


66C 


Ci 


283 


12. 


.758 


-528.936842 


-303.763297 


-273.223079 


-259.078765 


4. 


.83 


7.83 


66D 


c. 


274 


0. 


.640 




-301.866029 


-277.134660 


-273.454381 


7. 


.83 


9.74 


66E 


^ s 


274 


0. 


.573 




-301 464889 


-277.145223 


-273 573502 


9. 


.74 


11.50 


66F 


Cs 


274 


0. 


.552 




-301.193166 


-277.115736 


-273.602343 


11. 


.50 


24.31 


66G 


Cs 


271 


0. 


.033 




-299.437330 


-274.866840 


-271.669464 


24. 


.31 




67A 




318 






-544 754257 


-297 547886 










3.80 


67B 


Cs 


289 


15. 


.345 


-539.666842 [2] 


-308.979846 


-276.564507 




3. 


.80 


5.46 


67C 


C2 


288 


12. 


.945 


-538.179357[1] 


-309.130322 


-278.090470 


-263.853175 


5. 


.46 


7.72 


67D 




279 


0. 


.677 




-307 349356 


-282.171568 


-278 396059 


7. 


.72 


11.36 


67E 


C2v 


279 


0. 


.674 




-307.234099 


-282.166219 


-278.400953 


11. 


.36 


21.85 


67F 




276 


0. 


.033 




-304.922619 


-279.938518 


-276.682086 


21. 


.85 




68A 


^ s 


324 






-555 582496 












3.68 


68B 


Ci 


294 


15. 


.980 


-551.451174 


-314.374880 


-281.011142 


-265.753769 


3. 


.68 


7.29 


68C 




293 


12. 


.938 


-549.985909[1] 


-314.146641 


-281.905282 


-267.100885 


7. 


.29 


7.78 


68D 


Ci 


283 


0. 


.666 




-311.913082 


-286.249826 


-282.429781 


7. 


.78 


9.67 


68E 




283 


0. 


.597 




-311.515987 


-286.263045 


-282.554370 


9. 


.67 


11.24 


68F 


Cs 


283 


0. 


.573 




-311.246697 


-286.237292 


-282.588836 


11. 


.24 


11.51 


68G 




283 


0. 


.552 




-310 973538 


-286 207609 


-282.617642 


11. 


.51 


13.46 


68H 


C^y 


282 


0. 


.032 




-311.238864 


-285.954388 


-282.683003 


13. 


.46 




69A 


Ci 


329 






-566.364140 


-307.245602 










3.65 


69B 


^ 1 


308 






-560 639499 


-318 205833 






3. 


.65 


4.22 


69C 




298 


14. 


.909 


-559 783221 fll 


-319 687802 


-286.265131 




4. 


.22 


5.22 


69D 


c, 

^ s 


299 


16. 


.223 




-319.819905 


-285 836885 fll 




5. 


.22 


7.64 


69E 


c. 


288 


0. 


.705 




-317 396150 


-291.284253 


-287 366805 


7. 


.64 


9.66 


69F 


Ci 


288 


0. 


.653 




-317.293147 


-291.294050 


-287.462110 


9. 


.66 


18.06 


69G 


Ci 


286 


0. 


.034 




-315 861945 


-290 031196 


-286 695880 


18. 


.06 




70A 




335 






-577 739782 












3.22 


70B 


Ci 

i 


332 






-577 286914 


-313 362561 






3. 


.22 


3.67 


70C 




314 






-571.815209 


-324.211871 






3. 


.67 


4.16 


70D 


C^y 


304 


15. 


.324 


-570.747100 


-325.887749 


-291.872039 




4. 


.16 


7.89 


70E 


Cs 


293 


0. 


.676 




-323.082118 


-296.412149 


-292.439398 


7. 


.89 


12.18 


70F 


Ci 


293 


0. 


.659 




-322.802226 


-296.365659 


-292.462856 


12. 


.18 


17.96 


70G 




291 


0. 


.035 




-321 499078 


-295.108247 


-291 708765 


17. 


.96 




71A 


Ci 


330 






-588.396687 


-316.491227 










3.42 


71B 




320 






-583.202864 


-330.363241 






3. 


.42 


4.28 


71C 




310 


17. 


.608 


-581 944905 fll 


-331.588748 


-295 611936[2l 




4. 


.28 


7.29 


71D 




299 


0. 


.695 




-329.356130 


-302.405357 


-298.405353 


7. 


.29 


21.50 


71E 


Cs 


296 


0. 


.034 




-326.944961 


-300.178953 


-296.721356 


21. 


.50 




72A 


Ci 


335 






-599.690310 


-320.666142 










3.40 


72B 




314 


17. 


.960 


-595.064909 


-335.877119 


-299.300783 


-283.195227 


3. 


.40 


4.77 


72C 




313 


16. 


.818 


-592.641021 


-336.121753 


-299.750026 


-283.664026 


4. 


.77 


7.29 


72D 


Ci 


303 


0. 


.713 




-333.967611 


-306.460206 


-302.387256 


7. 


.29 


10.53 


72E 


Ci 


303 


0. 


.695 




-333.804475 


-306.453999 


-302.413229 


10. 


.53 


21.49 


72F 


c. 


300 


0. 


.038 




-331.710111 


-304.261423 


-300.734505 


21. 


.49 
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TABLE I. continued. 





PC 


n 


-^strain 


po — 3 


oi\ — 6 


On = 10 


na = 14 


Pmin 


Pniax 


73A 




339 






-610.936684 


-324.743483 










3.50 


73B 




319 


18, 


.316 


-603 873759 f2l 


-341.266253 


-304 026668 


-287 806170 


3 


.50 


6.97 


73C 


Cs 


318 


16 


.192 
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